NASA TECHNICAL NOTE 


04 

to 

CO 


i 




1 NASA TN H-^91 

C. ! 


LOAN COPY: F 


AFWL £ 


KIRTLAND AF 


ru : 
x: | 

IT • 


on 



EFFECTS OF TWO VIBRATION-DISSOCIATION 
COUPLING MODELS ON NONEQUILIBRIUM 
HYPERSONIC FLOW OF NITROGEN OVER A WEDGE 


by Fred R. Dejarnette 

Langley Research Center 
Langley Station, Hampton, Va. 



. rO • 


J < 
< 




WASHINGTON, 0. C. 






• AUGUST 1966 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECH LIBRARY KAFB, NM 



TECH library kafb, nm 

■Nunn 

013054b 

NASA TN D-3521 


EFFECTS OF TWO VIBRATION- DISSOCIATION COUPLING MODELS 

ON NONEQUILIBRIUM HYPERSONIC FLOW 

OF NITROGEN OVER A WEDGE 

By Fred R. DeJarnette 

Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — Price $2.00 



EFFECTS OF TWO VIBRATION- DISSOCIATION COUPLING MODELS 
ON NONEQUILIBRIUM HYPERSONIC FLOW 
OF NITROGEN OVER A WEDGE 

By Fred R. DeJarnette 
Langley Research Center 

SUMMARY 

The effects of the coupled vibration-dissociation (CVD) and coupled vibration- 
dissociation- vibration (CVDV) models on the inviscid, nonequilibrium flow over a wedge 
have been compared with the flow fields with no coupling (UVD model). The numerical 
method used for the flow-field computations makes use of Lax's finite- difference method 
but employs a modified- difference scheme on the wedge surface. Distributions of flow 
properties along and normal to the wedge surface are presented for the flow of nitrogen 
for free-stream Mach numbers of 15.25 and 23.76 past wedges with attached shock waves 
and wedge half-angles close to 43°. 

At the lower Mach number of 15.25, where dissociation is small, the CVD and 
CVDV models gave identical results. However, the UVD model predicted a much larger 
degree of dissociation near the tip of the wedge than the CVD and CVDV models. 

At the higher Mach number of 23.76, where dissociation is appreciable, large dif- 
ferences were observed between all three coupling models. The CVD model was found 
to have a dissociation overshoot in addition to the vibrational overshoot, whereas the 
CVDV model eliminated both of these overshoots. It was further found that the thickness 
of the shock layer is largest for the CVDV model and smallest for the UVD model. 

As a result of the present investigation, it is concluded that the vibration- 
dissociation coupling does have a significant effect on the flow fields past finite bodies 
at the higher Mach numbers. Even for the Mach numbers where dissociation is small, 
the effect of the coupling mechanism on the degree of dissociation is appreciable. 

INTRODUCTION 

A considerable effort is currently being devoted to the determination of flow fields 
surrounding reentry vehicles traveling at velocities and altitudes where the gas in the 
flow field deviates significantly from a perfect gas and where nonequilibrium real- gas 
effects must be considered. In addition to the usual nonequilibrium vibrations and 



dissociation associated with diatomic molecules, the coupling mechanism between the 
vibrational excitation and molecular dissociation has been shown to have a significant 
effect on the flow properties (refs. 1 and 2). Several coupling models have been devel- 
oped (refs. 1 to 5) and used in the computation of nonequilibrium flow fields (refs. 1 
to 13). The present paper examines the effect of two coupling models (compared with 
the uncoupled model) on the steady, inviscid flow of nitrogen past wedges. 

Hammerling, Teare, and Kivel (ref. 1) developed a coupled vibration- dissociation 
(CVD) model which accounts for the vibrational nonequilibrium effects on the rate of 
dissociation. However, the normal shock-wave flow fields computed in references 6 
and 7 showed that, at the higher shock strengths, the vibrational temperature overshoots 
the translational temperature. This phenomenon has been questioned on physical grounds 
in reference 14 and is no longer thought to be realistic. The main objection to the CVD 
model is that it does not account for the drain on vibrational energy due to the dissocia- 
tion process. This CVD coupling was modified by Treanor and Marrone (ref. 2), and 
independently by Heims (ref. 3), to include the converse effects of dissociation on vibra- 
tional relaxation. The coupled vibration-dissociation- vibration (CVDV) model eliminates 
the vibrational overshoot observed in the CVD model. 

Both the CVD and CVDV models contain two important assumptions: 

(1) The relaxation of the vibrational degree of freedom to the final Boltzmann dis- 
tribution at the equilibrium temperature is assumed to occur through a series of 
Boltzmann distributions. 

(2) Dissociation can occur with equal probability from any vibrational level in a 
sufficiently energetic collision. 

In reference 4 Marrone and Treanor modified assumption (2) so that dissociation 
could occur preferentially from the higher vibrational levels, and Marrone (ref. 11) com- 
pared the preferential model with the CVDV model for both normal and bow shock waves. 
A still more elaborate method was recently introduced by Treanor (ref. 5) for coupling 
the vibrations and dissociation by using non-Boltzmann distributions and different vibra- 
tion and dissociation rates for the individual vibrational levels. This model also accounts 
for the coupling between recombination and vibrational relaxation, which is important in 
expanding flows. Although these two latter modifications are more realistic, their 
application to flow fields (other than one- dimensional problems) is quite tedious. 

The purpose of the present paper is to compare the CVD and CVDV models with 
the uncoupled model (UVD) for the nonequilibrium flow of nitrogen past thick wedges with 
an attached shock wave. Although the CVD and CVDV models may oversimplify the actual 
coupling mechanism, they should be useful for assessing the importance of the nonequi- 
librium coupling effects on finite-body flow fields. 
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SYMBOLS 


M 

D 

E(T,T V ) 

Ey 

Ey,e 

f 

G(T) 

H 
H t 
• h 

K e 

k f 

k r 

l 

M 

m j 

N 

n 

P 


concentration of jth species, moles per cc 
dissociation energy per molecule 

average value of energy lost from vibration due to dissociation (per unit mass) 

vibrational energy per unit mass 

equilibrium vibrational energy per unit mass 

vibration- dissociation coupling factor defined by equation (26) 

average energy gained by vibration due to recombination (per unit mass) 

static enthalpy per unit mass 

total enthalpy per unit mass 

Planck's constant 

equilibrium constant for vibrational equilibrium 
specific forward reaction- rate constant 
specific reverse reaction- rate constant 
length scale 

locally frozen Mach number 
mass of a molecule of jth species 
concentration of nitrogen, moles per cc 

number of vibrational energy levels included in dissociation energy 
pressure, atm (unless otherwise noted) 


Q 


partition function 


R universal gas constant 

T translational temperature, °K 

T m parameter defined by equation (27), °K 

T v vibrational temperature, °K 

t time, sec 

u,v velocity components in x and y directions, respectively 

V velocity magnitude 

Wj molecular weight of jth species 

x,y coordinates along and normal to wedge surface, respectively 

a atom mass fraction, Pp^/P 

y locally frozen ratio of specific heats 

r\/y 

r rate of dissociation, 

9 wedge half- angle 

©D’®R’®v characteristic dissociational, rotational, and vibrational temperatures, 
respectively, °K 

dE v 

A vibrational energy rate, — — 

dt 

v vibrational frequency 

p mass density 

Pp dissociation density function given by equation (23) 
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r vibrational relaxation time, sec 

\p constant defined by equation (10) 

Subscripts: 

e equilibrium conditions 

el electronic 

f forward or dissociation reaction 

j species 

for general case: A or A 2 

for specific case of nitrogen: N or N 2 

r recombination or reverse reaction 

vib eq vibrational equilibrium 

°° undisturbed free- stream conditions 

Barred quantities are dimensionless (see eq. (7)); unbarred quantities are 
dimensional. 

Abbreviations: 

CVD coupling effect of vibration on dissociation 

CVDV coupling effects of vibration on dissociation and dissociation on vibration 
UVD uncoupled vibration and dissociation 

ANALYSIS 

The problem studied in this investigation is the effect of the CVD and CVDV 
coupling models on the steady, inviscid, hypersonic flow of nitrogen past thick wedges 
with an attached shock wave. A diatomic gas subject to harmonic vibrations, dissocia- 
tion and recombination is considered. Translation and molecular rotation modes are 
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assumed to be fully excited, whereas the electronic internal energy and ionization are 
neglected. For more detail concerning the gas model, see references 1, 2, 3, and 15. 

The numerical method used for the flow-field solution replaces the governing 
partial- differential equations with finite- difference equations. (See refs. 15, 16, and 17.) 
Then, with known flow properties along an initial data line near the tip of the wedge, the 
downstream flow field is computed by marching forward in the streamwise direction. 
This numerical technique is restricted to hyperbolic equations; hence, the steady- state 
flow field being computed must remain supersonic. 


Basic Equations 

The flow geometry and coordinate system are illustrated in the following sketch 
for a wedge. For this body-oriented coordinate system, x is measured along the sur- 
face of the wedge, y is normal to the sur- 
face, and the velocity components in these 
respective directions are u and v. The 
basic flow and rate equations (specialized 
for a wedge) are written in conservation 
(or divergence) form as (ref. 15) — 


y.v 



y-momentum: 


Continuity: 


8 (pu) , 8 (pv) _ Q 

dx 9y 


x-momentum: 


8 (p + pu 2 ) 8(puv) _ Q 

9y 


dx 


8 (pu v) a(p + pv 2 ) _ 0 
dx dy 


( 1 ) 


( 2 ) 


(3) 


Vibrational energy rate: 


~a(p U Ey) 9 (pv Ey) 

Qy 


p* w a 2 a 


9x 


V^R 


= 0 


(4) 


Rate of dissociation: 


a(pua) 9(pvof) £_IT _ n 
9x 9y " 


(5) 
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Energy: 


—2 —2 

H + u * v = H t * Constant 
2 


( 6 ) 


The non dimensional vibrational energy is Ey, a is the atom mass fraction (p^/p)? 

A is the dimensional vibration rate (dE v /dtj, and r is the dimensional dissociation 
rate (do/dt). All quantities have been nondimensionalized as follows (barred quantities 
are nondimensional): 


u.v 

u, v = — 2 — 

* V 

v oo 


p = p 


T = 


E v = 


H = 


PooVoo" 

H 


D = - e _ 


x,y = ^ 


e v wa 2 


RT„ 


s s s _ 0 D’°R’ 0 v 
0 D’ 0 R>°v “ ™ 


The equations for enthalpy and state (see refs. 15 and 16) are — 
Enthalpy: 


H = 


"j^7_H-_3 


+ (1 - O') E v + ctO 


D ] 


State: 


where 


P = pT(l + a)ip 


* s 


RT 0 


w a 2 v » 2 


= [r»(i + oo)m„ 2 ] 


The locally frozen ratio of specific heats is 

7 + 3o! 


r = 


5 + 0! 


and the locally frozen Mach number is 


M = 


p ♦ vA 1/2 

(u 2 + v 2 )wa 2 

\ rP/p / 

y( 1 + q;)RT 


1/2 


(7) 


( 8 ) 

(9) 

( 10 ) 

(ID 


( 12 ) 


The coupling models affect the vibrational and dissociational rates, A and r. The 
characteristic dissociation temperature for nitrogen, 0jj, is 113,240° K. 
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Vibrational Rate 


A linear harmonic oscillator model is assumed for the diatomic vibrations. The 
vibrational relaxation, neglecting the dissociation coupling effect, is from reference 18 

B A = E v,e ~ E v (13 

dt r 

where the equilibrium vibrational energy at the local translational temperature is 


W A2 ie 0 v/ T - ]) 

The characteristic vibrational temperature 0 V for nitrogen is given by reference 19 
as 3396° K. Reference 2, where numerical fits were made to Blackman's data (ref. 20), 
gives the vibrational relaxation time r in sec for nitrogen as 


r = (l. 


1 x 10 ' 


t 


Equation (13) represents the vibrational energy rate for both the UVD and CVD 
models. However, it must be modified by an expression which accounts for the coupling 
effect of dissociation on vibrational relaxation when considering the CVDV model. 

Dissociation Rate 

For the dissociation of nitrogen, reference 2 gives the following chemical 
reactions: 

*1 

N 2 + N~2N + N (16) 


N 2 + N 0 — »- 2N + N 2 

kr 2 

For nitrogen, the experimental forward reaction- rate constants are not known; 
however, reference 2 gives the reverse reaction- rate constants as 

k r = (3.33)10 21 T' L5 

mole -sec 


k r „ = (3.33)10 19 T -1 - 5 

mole -sec 
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The net rate of production of the atom mass fraction for the reactions given by equa- 
tions (16) and (17) is (ref. 15) 


2kj*jPQ! 

r*»i (1 ■ “ ,wn * . j 

k r2 p(l - a) 

kf 2 (1 - q)W N2 

Wn 2 

k ri 4 P 

W N 2 

k r2 4 



( 20 ) 


Equation (20) gives the dissociation rate once appropriate expressions for the forward 
reaction-rate constants have been determined. These constants are in turn dependent 
on the particular coupling considered (UVD, CVD, or CVDV). 

When the vibrational energy is in equilibrium with the translational mode, the for- 
ward and reverse reaction- rate constants are related through the equilibrium constant 

Re = (21) 


(It is assumed that vibrational nonequilibrium does not affect k r .) The equilibrium 
constant K e can be calculated from molecular properties by the use of quantum statis- 
tical mechanics. For the reactions given by equations (16) and (17), the equilibrium 
constant is (ref. 21) 


K e = 


4 P D e 


-0 D /T 


Wn 2 


( 22 ) 


where 


m N 


P D " 


m N( 7rm N kT ) 
2h 3 


3/2 





mass of a molecule of atomic nitrogen 


(23) 


h Planck's constant 


k Boltzmann’s constant 

characteristic rotational temperature, 2.847° K for nitrogen 
Q e j electronic partition function 
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As shown in appendix B of reference 19, the electronic partition functions for 
nitrogen vary less than 5 percent for 1500° K < T < 8000° K. Therefore, average 
values over this temperature range are just the degeneracies of the ground states 

Q eZ,N 2 = 1 Q eZ,N = 4 


Uncoupled Vibration- Dissociation (UVD) Model 

In the earlier analyses of nonequilibrium flows the vibrational relaxation and 
molecular dissociation were not coupled. The measured reaction- rate constant was 
used as the actual reaction-rate constant, and the opposing reaction- rate constant was 
determined through the equilibrium constant. Thus the UVD model sets 

kf ~ (kf) v ib e q 


and the forward reaction- rate constants for the UVD model are (from eq. (21)): 


kfj - K e k ri 

(24) 

k f 2 = K e k r 2 

(25) 


The vibrational rate equation for the UVD model is given by equation (13). 


Coupled Vibration- Dissociation (CVD) Model 

The use of the equilibrium constant to relate the forward and reverse reaction- 
rate constants (eq. (21)) is justified only when vibrational equilibrium exists. 
Hammerling, Teare, and Kivel (ref. 1) developed the CVD coupling to account for the 
effect of vibrational relaxation on molecular dissociation. By using the harmonic- 
oscillator calculations of Montroll and Shuler (ref. 22), the ratio of the dissociation- rate 
constant to the value with vibrational equilibrium was found to be 


where 


k f _ l| 

\ _ e - n0 v/ T m\ 

(e 0 v/ T v - l) 

£ 

cr 

CD 

►Q 

1 

PI 

V e 0v / Tm - 1 / 

(e°v/ T - l) 


1 _ 1 1 _ 

T m = T v " T 

R0y 

Ev = ^Vt7^J 


(26) 


(27) 

(28) 
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and n is the smallest integer greater than D/ht>, that is, integral number of vibra- 
tional levels up to the dissociation cut-off, which for nitrogen is 34. 

Thus the nitrogen reverse- rate constants for the CVD model are given by equa- 
tions (18) and (19), whereas the forward reaction- rate constants are 

kf l = f ( kf l)vib eq = 1 kr l(M (29) 


kI 2' f N) vibeq = £ MM 


(30) 


The effect of the coupling given by equation (26) is to slow down the dissociation 
rates in the region where E v < E v e since 0 < f < 1. However, in the region of the 
vibrational overshoot (t v > T), which occurs with the CVD model, the dissociational 
rates are increased because f > 1 for T v > T. 


Coupled Vibration- Dissociation- Vibration (CVDV) Model 

The CVD model does not include the converse effect of molecular dissociation on 
vibrational relaxation. To account for this additional effect in the coupling of vibrations 
and dissociation, Treanor and Mar rone (ref. 2) modified the vibrational- energy- rate 
equation given by equation (13) to the form 



where is the concentration (moles/cc) of the diatomic molecules considered, and 

the subscripts f and r refer to the forward (dissociative) process and reverse (recom- 
bination) process. The e(T,T v ) term is the average value of the energy lost from vibra- 
tion due to dissociation, G(T) is the average energy gained by vibration due to recom- 
bination, and the last term in equation (31) accounts for the change in the number of oscil- 


lators during dissociation. Note that both 
quantities in equation (31). 



are used as positive 


Using the same assumptions as for the CVD model, Treanor and Marrone obtained 


R0 V 

1 

n 

w N2 

( e 0v/ T m . 

( e n0 v/ T m _ i) 


(32) 


11 



and 



Therefore, the vibrational rate for the CVDV model is 




( 33 ) 



The same reaction- rate constants are used for the CVDV model as for the CVD model. 
For the reactions given by equations (16) and (17), 



At relatively low shock strengths, where the degree of dissociation is small, the 
CVD and CVDV models give identical results. At higher shock strengths, the CVDV model 
eliminates the vibrational temperature overshoot observed in the CVD model. The CVDV 
model in turn predicts a slower rate of dissociation since T v < T (and thus f < 1) 
throughout the entire relaxation process. 

Method for Numerical Computations 

The numerical method used for the flow-field computations is described in refer- 
ence 16. Basically, the method replaces the partial- differential equations with finite- 
difference equations, by using the difference scheme of Lax (ref. 17) along with a 
modified- difference method for the boundary points. As a consequence of Lax’s difference 
scheme an artificial viscosity is implicitly introduced and allows the computations to pro- 
ceed downstream of an initial data line as if no shock wave were present at all. The shock 
wave appears in the solution, however, smeared over several mesh spaces while accu- 
rately giving the proper jump conditions across the shock. Approximate flow properties 
on the initial data line, which is near the tip of the wedge and normal to the surface, are 
determined from the wedge-tip gradients given in appendix C of reference 16. 
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Since the system of equations being solved is hyperbolic, the ratio of the mesh 
spacing in the x-direction to the mesh spacing in the y-direction is limited by the hyper- 
bolic stability criterion. The free stream as well as the shock layer must be considered 
in determining the stability criterion because the computations cariy through the shock 
wave. For all the examples given here, the free- stream stability criterion was more 
stringent than that of the shock layer since the free-stream velocity is inclined at the 
wedge half -angle 9 to the coordinate system. (See ref. 16 for the details.) 

RESULTS AND DISCUSSION 


References 1, 2, and 12 show the effect of the CVD and CVDV couplings on the one- 
dimensional flow downstream of normal shock waves. However, these effects have not 
been investigated for flow fields past bodies. Therefore, the following cases were com- 
puted to determine the differences between the UVD, CVD, and CVDV gas models for flow 
over wedges. In each case the shock-wave angle at the tip of the wedge was chosen to be 
59°. The Mach numbers were chosen so that the shock strengths at the wedge tip would 
be the same as those for the normal shock solutions computed in reference 2. 

Case I 


Gas = Nitrogen = 15.24987 6 = 43.066° 

Too = 285° K P^lmmHg E voo = 0 

Oioo = 0 

Case n 


Moo = 23.75726 

Poo = 1 mm Hg 

Figures 1 to 4 give the variation of P, T, E v , and a along the wedge surface 
for Case I (Moo = 15.24987). These results follow the trend shown in reference 16; and, 
to the scale of these figures, there is no difference between the CVD and CVDV models. 
The UVD model gives the same results as the CVD and CVDV models for P, T, and 
Ey However, figure 4 shows that the degree of dissociation is much larger for the UVD 
model than for the other two, particularly near the tip of the wedge. Due to the small 
magnitudes of a, these differences have very little effect on the other flow variables. 


Gas = Nitrogen 
Too = 285° K 

Ooo = 0 


6 = 43.319° 
E v ,oo = 0 


Figure 5 gives the pressure distribution along the wedge surface for Case H 
(Moo = 23.75726). In order to carry the solutions out to x = 100 cm, it was necessary to 
increase the mesh scale for x > 2 cm. After increasing the mesh scale, the behavior of 
the pressure was questionable, but still maintained a value close to that at x = 2 cm. 
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Therefore, the pressure curve in figure 5 was stopped at x = 2 cm. Although the pres- 
sure was questionable for x > 2 cm, the other properties should be accurate. The justi- 
fication for this fact is that the comparisons given in references 16, 23, and 24 show that 
accurate surface properties can be obtained even where the pressure is not accurate. 

Figures 6, 7, and 8 illustrate the large differences between the UVD, CVD, and 
CVDV models for T, E v , and ou In addition to the vibrational overshoot for the CVD 
model (which can be deduced from fig. 7), figure 8 shows a dissociation overshoot 
(a; > cte) for the CVD model. This dissociation overshoot has not previously been 
reported in other flow-field computations. As a consequence of the vibrational and dis- 
sociational overshoots, the T and E v curves for the CVD model exhibit "buckets" in 
their surface distribution. The CVDV model eliminates both of these overshoots. Fig- 
ure 8 shows the same early increase in a for the UVD model as was found in Case I. 

The profiles of the flow properties normal to the wedge surface at x = 2 cm are 
given in figures 9 to 12 for Case II. All these figures show that the shock layer is largest 
for the CVDV model and smallest for the UVD model. 

The large effects produced by the coupling models in Case II can be explained by 
the coupling factor f given in equation (26). A study of the results obtained here 
revealed that the forward reaction rate dominated the reaction over the major portion of 
the flow field downstream of the tip of the wedge. The forward- rate constant was deter- 
mined from 

k f = f(K e )k r 

and for the UVD model f b 1, whereas for the CVD and CVDV models 

0 < f < 1 for E v < E v> e 
f > 1 for E v > E v>e 

This causes the CVDV model dissociation rate to be slower than for the UVD model since 
E v < Ey ?e throughout the flow field. However, the CVD model has E v < E v e near the 
tip of the wedge, but the vibrational overshoot (E v > E v e j makes f > 1 which in turn 
increases the dissociational rate over that predicted by the UVD model and leads to the 
dissociation overshoot observed for the CVD model. 

It should be noted that, for the temperatures and pressures encountered in Case II, 
the effects of vibrational anharmonicity, electronic internal energy, and ionization are 
not necessarily small, although the present analysis does not consider these effects. 
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CONCLUDING REMARKS 


The effects of the CVD and CVDV coupling models on the nonequilibrium flow of 
nitrogen past wedges have been compared with the uncoupled (UVD) model. Mach num- 
bers of 15.25 and 23.76 were considered for wedges with half-angles near 43°. 

At the lower Mach number of 15.25, where dissociation is small, the CVD and 
CVDV models gave identical results. However, the UVD model predicted a much larger 
degree of dissociation near the tip of the wedge than the CVD and CVDV models. 

For the higher Mach number of 23.76, where dissociation is appreciable, large dif- 
ferences were observed between all three coupling models. The CVD model was found 
to have a dissociation overshoot in addition to the vibrational overshoot, whereas the 
CVDV model eliminated both of these overshoots. It was further found that the thickness 
of the shock layer is largest for the CVDV model and smallest for the UVD model. 

As a result of the present investigation, it is concluded that the vibration- 
dissociation coupling does have a significant effect on the flow fields past finite bodies at 
the higher Mach numbers. Even for the Mach numbers where dissociation is small, the 
effect of the coupling mechanism on the degree of dissociation is appreciable. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., April 6, 1966. 
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Figure 2.- Temperature along wedge surface for Case l. 
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